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Detecting gravitational wave bursts (characterised by short durations and poorly modelled wave- 
forms) requires to have coincidences between several interferometric detectors in order to reject non 
stationary noise events. As the wave amplitude seen in a detector depends on its location with 
respect to the source direction and as the signal to noise ratio of these bursts are expected to be 
low, coincidences between antennas may not be so likely. This paper investigates this question from 
a statistical point of view by using a simple model of a network of detectors; it also estimates the 
timing precision of a detection in an interferometer which is an important issue for the reconstruction 
of the source location, based on time delays. 
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I. INTRODUCTION 

In the next few years, the first generation of long baseline interferometric detectors of gravitational waves (GW) 
will be operational. Among the most promising sources of GW, the compact binary inspirals and the periodic 



OC ■ sources have been already studied for a long time (see [^ or ||6| for a review); more recently, some analysis methods 
Li^ ,. have been developed to look for GW bursts of short duration and unknown waveform which are the subject of this 
O ■ paper. 



^yi\ Type II supernovae ( [Q and references therein) and the merging phase of binary compact star and/or black hole 



systems M are the most common burst sources foreseen but others such as e.g. cosmic strings p| could also be 
•^ ■ considered promising. The lack of knowledge of such GW signals prevents from using the optimal (Wiener) filtering 
5^ method and requires sub-optimal (i.e. less efficient yet robust) filters; various techniques have already been studying 



by different groups. The "excess power statistics" |10 monitors the power of the detector output along the time and 
is shown to be optimal when minimal hypothesis (signal duration and bandwidth are the only known quantities) are 
taken into account [Q. Time- frequency methods have also been presented: the authors in [Q use the Wigner-Wille 
distribution to transform time series into two-dimensional maps in which ridges are looked for by using the Steger 
algorithm; a Student-like test on noise periodograms has been studied in order to detect non-stationary events [|l3| . 
More generally, we have developed a set of filters able to be used as on-line triggers [[l4|-p^ , the most promising being 
the "alternative linear fit filter" (ALF) |T^, which aims at detecting a nonzero slope in data windows. 

As burst signals are poorly modelled, only an accurate knowledge of the detector behaviour could allow to properly 
separate a real GW from non-stationarities in a single antenna; therefore, one will be more confident about the reality 
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of a candidate event in case of a significant coincidence between several interferometers, or between an interferometer 
and other types of detector (neutrinos, 7-rays satellites...). Many studies in the literature assume such multi-detections 
as the starting point of the analysis and deal with the "inverse problem" , i.e. how to extract informations - source 
location in the sky, GW waveforms, astrophysical parameters - from these triggered data |17|- p^ , pl[ . In this article 
we will not address this question which will be the main topic of forthcoming papers aiming at developing coincidence 
strategies for the detection of bursts with a network of both GW antennas and detectors sensitive to other radiations 
such as neutrinos. 

Simulations of massive star collapses [ P2|j23[ give typical GW amplitudes too low to expect a likely detection of 
sources outside the Galaxy; more optimistic estimates have been computed recently for the merging phase of two 
neutron stars [24) or a neutron star and a black hole [g5|, whose coalescences could be detected until 10 Mpc given 
the planned sensitivity of the first generation of detectors. In both cases, as the GW amplitude is just above the 
noise level, the signal to noise ratios (SNR) are nevertheless expected to be low, all the more that the wave amplitude 
detected in a given antenna depends on its location and orientation with respect to the source [p[: when the source 
direction is not orthogonal to the detector plane, the response to the GW is not optimal and so the SNR is decreased. 

Therefore the feasibility and the efficiency of coincidences between interferometers currently under construction 
worldwide is an important question: are these coincidences likely or not? This problem has attracted so far little 
attention in the literature whereas it should be the starting point of the study of coincidence strategies: a precise 
knowledge on how a network of detectors surveys the sky is essential in defining the best possible methods to analyse 
their data. 

This paper addresses this topic from a statistical point of view by performing Monte-Carlo simulations of the 
detection process in a network of antennas. After having described the model used in the computations, we briefly 
recall the structure of the beam pattern functions - describing the interaction between a GW and a detector - and 
extract from them some information on the accuracy with which a particular sky direction is surveyed by a given 
interferometer. We will use the celestial sphere coordinates, i.e. the right ascension a and the declination 5 to label 
sources. Sky maps will be constructed to characterise the detection power of each antenna and any combination of 
them. 

Next, we introduce a quantity independent of any detector, Pmax, characterising a GW strength; the simulation 
results - detection probabilities for different configurations in the network: single antenna, twofold coincidence, 
threefold coincidence - are then presented either as function of pmax and {a, S) on a sky map, or simply versus Pmax 



after averaging over all the sky directions. 

Finally, the timing accuracy is strongly connected to the detection problem: reconstructing the source location in 
the sky or performing coincidences with other types of detectors (7-rays, neutrinos...) requires a good knowledge of 
the GW event timing in the interferometer(s). Therefore, we also present some results about the timing precision one 
can expect for typical burst detection before concluding. 

II. MODELLING THE DETECTION PROCESS IN THE NETWORK OF INTERFEROMETERS 

In this section we present the different hypothesis used in the Monte-Carlo simulations whose results are presented 
below: first about the network of detectors; then we give some assumptions about the GW burst and the calibration 
of its amplitude; finally, we detail the detection process itself. 

A. The network 

Each detector Di is described by a set of four angles: its latitude li, its longitude Li, the angle between the two 
arms Xi ^nd 2111 orientation angle 7^ defined as the angle between the local South direction and the bisecting line of 
the detector arms counted counterclockwise. The values of these angles can be found in Table I extracted from [E6[ . 



Detector 


Latitude / 


Longitude L 


Arms 'separation' x 


'Azimuth' 7 


VIRGO 


43.6 


-10.5 


90.0 


206.5 


LIGO Hanford 


46.5 


119.4 


90.0 


261.8 


LIGO Livingston 


30.6 


90.8 


90.0 


333.0 


GEO600 


52.3 


-9.8 


94.3 


158.8 


TAMA300 


35.7 


-139.5 


90.0 


315.0 



Table I: detector data, all angles given in degrees 
The interferometer noise components are assumed to be white (after whitening |p^,E8| and applying line removal pfl] 
methods), Gaussian, broadband, additive and uncorrelated; lacking realistic sensitivities validated by measurements, 
- these informations will appear when detector runs start - we also assume they all have same noise standard deviation 
a ||30|] and sampling frequency /q. 

B. The gravitational wave burst and the source direction 



As a GW signal, we use an one-parameter set of Gaussian bursts 



G.(i) = ifexp(- — 



(2.1) 



where the width r is taken between 0.1 ms and 10 ms and X is a scahng factor. Therefore, the maximum SNR 
(corresponding to an optimal detector orientation with respect to the source direction p^ and by using Wiener 
filtering), referred as Pmax in the following, is given by Q: 

Pmax = J< (^-^j 



By assuming a value for pmax, it is straightforward to compute K fronr Eq. ( |2.2| ) and so to characterise the amplitude 
of the GW burst independently of any detector. 

In addition to the two coordinates a and 6 previously introduced, a third variable is necessary to determine 
completely the GW: the polarisation angle ip, one of the Euler angles describing the wave coordinate system in the 
TT gauge with respect to the celestial frame. In the simulations, unless specified otherwise, the GW sources are 
assumed to be uniformly distributed over the sky (a £ [— tt; tt], sin (5 G [—1; 1], ip G [— tt; tt]). 

Finally, one has to notice that the arrival time of the GW in the detectors is a priori not synchronised with the 
sampling which will cause some losses in SNR and in timing accuracy for very narrow bursts. 

C. The detection procedure 

The detection method used is Wiener filtering - correlation with the known signal itself - for sake of simplicity: 
it is an approximation of no consequence as an earlier paper mM showed that the one-dimensional parameter space 
[Tmin = 0.1 ms; Tmax = 10 ms] Can be covered by a discrete lattice of Gaussian filters ensuring a detection with a 
mismatch SNR loss lower than 1 %. The threshold value on the filter output is set to 77 = 4.89 which corresponds to a 
false alarm rate of 10~^, i.e. 72 per hour for the VIRGO value of the sampling frequency /q. One claims a detection if 
the filter output overcomes the threshold at least once in the analysis window. Decreasing the false alarm rate (10~^, 
10^^...) will not dramatically increase the threshold, so the results presented in the following remain essentially the 
same. 

III. AVERAGED BEAM PATTERN SKY MAP 
A. Beam pattern functions 

In the frame associated to the GW (wave direction being z by convention), in the TT gauge, the spatial metric 
perturbation is given by 

H{t)^ \ h^{t) -h+{t) (3.1) 

0/ 



with /i_|_ and h^ corresponding to the two independent wave polarisations. By assuming the detector size small 
compared to the reduced wavelength of the GW g and by following the treatment presented e.g. in jl^,^ it is 
straightforward to show that the response h{t) of an interferometric detector to this wave is a linear combination of 
the two polarisations: 

h{t) = F+{t)h+{t) + FAt)hy,{t) (3.2) 

The corresponding weighting factors are called beam pattern functions; they have values in the range [-1;!] depending 
on the longitude and the latitude of the detector location, as well as its orientation, the angle between the interferometer 
arms x, the sky coordinates (a, 5) of the source and the wave polarisation angle ?/'■ Still following the notations and 
the analysis of |p6| one can compute the general expression of the beam pattern functions: 

FAt)) ~ ^l^-sin2V' cos2V^j \b{t) ) ^'^■'^> 

The factor sin x recalls that the best response is achieved for detectors with orthogonal arms (like for instance VIRGO 
and the two LIGO interferometers); the a and h factors depend neither on V' nor on x and so 2?/' appears like a rotation 
angle and sinx like a scale factor for the beam pattern functions. 

Due to the Earth proper rotational motion the sky is in apparent motion with a period equal to a mean sidereal 
day. Therefore, the beam pattern functions associated to a source location (a, (5) depend also on the UT time i; let 
us introduce the local sidereal time T{t) for a detector of longitude L: 

T{t) = nt + TGrccnwich(O) - L (3.4) 

with K « 1.0027379 x 157hour 

rGiocnwich(O) is the Grccnwich sidereal time at Oh UT and the minus sign before L comes from the fact that longitudes 
are counted positive westwards. Then one can define the local hour angle of the source: 

f){t) = T{t) -a = Kt-{a + L)+ TcrccnwichlO) (3.5) 

Extensive calculations yield to compute the complete expressions of a and b (depending on a, d, I, L, 7 and i); we 
recall them for sake of completeness even if they are completely equivalent to those presented in section 2.1 of pfl] 
(see also references therein): 



a{t) = ^ T7 sin 27 (3 — cos2Z)(3 — cos 2(5) cos2Sj{t) — - cos 27 sin/ (3 — cos 25) sm2Sj{t) 
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— -sin 27 sin 2/ sin 2(5 cosSj{t) — -cos 27 cosl sin2(5sinio(t) — -cos 27 cos^ / cos^ S (3-6) 

b{t) = — cos27 sinZ sin (5 cos2i3(t) H — sin 27 (3 — cos2Z) sin (5 sin2i3(i) 

— cos27 cosl cos (5 cosi3(i) + - sin 27 sin 2/ cos5 smSj{t) (3.7) 



B. Detectors sky maps 

The main interest of studying beam pattern functions is to characterise the sensitivity of an interferometric detector, 
given for different directions in the sky, independently of its noise curve. This can be done by performing a quadratic 
average of the beam pattern amphtude over the unknown polarisation angle ■(/'• It leads to a f/^-independent quantity 
Fit) m given by 



Fit)^'-^^a^t)+b^t)e 



-T. 



(3.8) 



We choose to present the variations of i^ in a two-dimensional contour plot figure - the "sky map" - where any 
direction is located by a couple (a, sin 5) G [— tt; tt] x [—1; 1]. For a given detector, F depends not only on a and 5 but 
also on the time t. For sake of simplicity, each sky map shown below in the paper is determined at the time t = to 
defined by Kio + 2l3rconwich(0) — [2tt]; to have results at another time t one can simply imagine the map on a cylinder 
whose axis is along the sin (5 direction and to rotate it by SSj{t) = Sj{t) — Sj (to). 

Figure H compares the sky maps of the ?/;-averaged beam pattern function F for the five interferometers currently 
in development in the world: VIRGO, the two LIGO antennas, GEO600 and TAMA300. In each case, the F averaged 
value over the whole sky is 0.42 and the RMS is 0.16; each map presents 2 maxima corresponding to directions which 
are orthogonal to the detector plane and 4 minima (direction of the bisector of the two arms and its three images by 
rotations of angles 7r/2, n and 37r/2 in the detector plane) included in a large "valley" where F remains quite small. 

To simplify the analysis presented in this paper we now consider only three of these antennas as this is the minimum 
number required in order to be able to reconstruct the source direction in the sky in case of a full coincidence. We 
keep the two LIGO interferometers as they have been planned to work coherently by construction and as they have 
the same design sensitivities; VIRGO and GEO600 have almost identical F sky map (at least for the areas where 
the averaged beam pattern function is maximal) and so keeping one of the two (retaining VIRGO which is expected 
to have a sensitivity similar to LIGO) is enough to compute the main properties of an American-European network. 
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Finally, despite that TAMA300 sky map is different from the other ones, its expected sensitivity is too small - with 
respect to the other antennas previously mentioned - to help the network detection. So, in the following, the network 
chosen is made of three detectors: VIRGO (V) and the two LIGO interferometers in Hanford (H) and Livingston 
(L); the selected antennas are those expected to achieve the best broadband sensitivities. Nevertheless the results 
presented in this paper can be easily extended to accommodate additional interferometers into the network. 

The regions of maximum F are rather close for the two LIGO interferometers but quite different from those of 
VIRGO. Therefore one can expect on one hand a good complementarity between the three interferometers for a single 
detection - a clear advantage to maximise sky coverage. On the other hand, the ability to perform coincidences 
between the 3 detectors is expected to be significantly reduced - at least for small SNR values. 

C. Daily averaging 

As mentioned before, the two-dimensional maps previously defined are also function of time. Therefore, another 
way to present the information they contain is to average F over one sidereal day (and so over a) to keep only the 
dependence in the declination variable sin (5. 

(F) "™''"^''' (sin d) ^ ^ r daF{a, sin 6) (3.9) 

Figure presents the graphs of (_F) for the three detectors VIRGO, LIGO Hanford and LIGO Livingston. 

One can note that the vertical axis scale is zero-suppressed to enhance the variations which are in fact small (< 30%) 
around the mean value of 0.42. Due to the successive averages, these plots mainly depend on the detector's latitude 
which explains why VIRGO and LIGO Hanford curves are quite similar. To make this figure more concrete, vertical 
lines are given corresponding to the location of some centres of galaxies where GW sources are expected to be found: 
the Galaxy, the Magellanic clouds, M31 Andromeda and M87 in the Virgo cluster. For the galactic centre and M31, 
the three detectors are roughly identical whereas Livingston is better for M87 and VIRGO/Hanford for the Magellanic 
clouds. 

IV. DETECTOR CONFIGURATIONS 

The detection process (exceeding or not a threshold) is not linear; therefore, it can only be studied by using Monte- 
Carlo simulations. So, in the following, no a priori average on -0 is performed; we assume uniform distributions of 
sources in the sky which corresponds to the ranges of variables quoted in Tabic II. The results are twofold: 
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• Detection efficiencies versus pmax averaged over a, sin 6 and V'; 

• Detection sky maps, i.e. detection probability versus (a, sin (5) for a given Pmax, averaged on tp. 

Various configurations are studied, from a single interferometer to the search for full coincidences between VIRGO 
and the two LIGO antennas. 

A. Single detector 

1. The example of VIRGO 

As shown in section III-A, the amplitude of the GW interacting with a detector has to be multiplied by the beam 
pattern function with respect to the ideal case (optimal incidence and polarisation). Therefore the signal to noise 
ratio is lowered by the same factor: let us call prcai the product F x pmax where F is the beam pattern function. 
Figure presents the detection efficiency as function of the real signal to noise ratio Prcai for a linear filter with a 
10^^ false alarm rate |p3. The 50% efficiency is reached for pioai = threshold w 5 and the detection is almost always 
successful for real SNR higher than 7 (efficiency higher than 98%). Note that this curve obviously does not depend 
on the beam pattern. 

Adding the effect of the non optimal detector response gives graphs such as those shown on Figure 0. For two 
values of Pmax - 10 and 20 - and a Gaussian width of t = 1 ms, one sees 

• continuous line: the distribution of Pieai for Pmax constant; this curve shape is only due to the beam pattern; 

• dashed line: the fraction of events really detected, computed by including detector efficiency from Figure 0. 

From to £i^ the distribution of pioai is flat; for higher values, it decreases monotonously up to the maximal signal 
to noise ratio. Adding the noisy detection process considerably lowers the number of events detected at small real 
SNR but keeps almost all those with pioai > 7. By comparing the number of detected events with the total number 
of events generated by the Monte-Carlo simulation, one can compute the detection efficiencies: 32% for pmax = 10, 

66% for Pmax = 20. 

Simulating the real detection process, one obtains Figure ^ which shows the detection probability versus Pmax for 
different values of t. As expected, this probability only depends on the signal to noise ratio which is the relevant 
variable in signal analysis. Due to the beam pattern functions, the detection efficiency remains low even for large 
optimal signal to noise ratios: with Pmax=10 the detection efficiency is only about 30%, it reaches 50% for Pmax=14 
and 90% for - unlikely ~ high values. 



From those graphs one can conclude that including beam pattern functions in data analysis has major consequences 
on detection efficiency by reducing greatly the detection probability even for high Pmax- 

The previous results deal with averaged efficiencies around the whole sky; of course, the detection probability 
depends on the relative position of the source with respect to the interferometer. Figure ra presents the detection 
efficiencies versus (a, sin (5) for three values of /Omax- One can note that these three maps are highly correlated with 
the upper graph of Figure ^: F is thus a good - and easily computable - estimator of the detection probability in a 
given sky direction, especially in the low SNR region where are likely to be located the first detected events. 

In the ffi'st case (pmax = 10) detections are possible for incoming wave directions almost orthogonal to the detector 
plane reducing considerably the detection rate: a large fraction of the sky remains invisible. For pmax = 15 and 
Pmax — 20, the detection probability is more uniform over the sky; the four minima are still there but they cover 
now narrower regions whose areas decrease as Pmax increases. In the last case, there are simply blind islands among 
regions where the detection efficiency is larger than 60%. 

Converting the results of these maps into fractions of sky with given detection probability leads to Figure n One 
can note that especially at large Pmax the fraction of sky covered decreases first slowly when the detection probability 
level increases and then more rapidly after some corner value. The reason for this evolution can be understood by 
looking at Figure pi considering the large part of the sky where the efficiency is large and almost constant. The 
detection probability is higher than 30% in 40% of the sky for Pmax=10; for Pmax=15 it is more than 50% in half of 
the sky and it is almost 60% in 70% of the sky for pmax^SO. Therefore detections are likely for high maximal SNR 
Pmax in any direction but even for pinax=20, the efficiency never reaches 90%. 

2. Daily averaged detection probability 

This analysis can also be done for the two LIGO interferometers; as expected, averaging the detection probability 
over the whole sky while taking pmax constant gives exactly the same results. To compare the three interferometers, 
we simply average on a over one day. So, Figure H shows the daily averaged detection efficiencies for each of the 
detectors and for three different values of the maximal signal to noise ratio, pmax ~ 10, 15 and 20. As on Figure 0, the 
dashed vertical lines indicate some locations of galaxy centres. One can note that the detection probability is shifted 
to higher values for increasing Pmax! nevertheless, the curve shapes do not change too much and remain close to those 
shown on Figure for the averaged beam pattern functions: LIGO Livingston has the best results for small values 
of I sm6\ whereas LIGO Hanford and VIRGO are more efficient for large | sin(5|. But there remain some differences 



between the two Figures so that one should refine our previous statement: computing averaged beam pattern maps 
allow to have an idea of the detection probability in different sky areas but it is necessary to perform Monte-Carlo 
simulations including the detection process to compute the correct probabilities. 

B. Detector complementarity 

As we already checked that detection depends only on signal to noise ratio and not on the signal shape itself, we 
choose T—1 ms - a typical burst duration for the Gaussian waveform we use - in the following. In this section, we 
study the complementarity of the antennas, i.e. how their sky coverages complement each other, cither for a single 
detection - the "OR" strategy - or for different types of coincidences - the "AND" strategies. 

Figure ^ shows the detection probability versus Pmax for different configurations of detectors in the network; the 
curve already presented on Figure g| - single detector efficiency - is simply recalled for comparison. The continuous 
line shows the efhciency of detection in any single detector, i.e. the probability for a detection in at least one of the 
three interferometers. Network detection potential is clearly better than for a single antenna and the smaller the /Omax 
the higher the difference: more than twice more detections for 5 < Pmax l£ 10, more than 1.5 below Pma.K=^'^- The 
50% efficiency is reached at Pmax^S and for pmax > 15 the probability is higher than 85% on average, corresponding 
to a likely detection in most parts of the sky. 

On the same Figure, one can also see the curve corresponding to a coincidence detection in at least two detectors 
among three ("OR" of twofold coincidences). For small values of pmax, the twofold coincidence probability remains 
lower than the detection efficiency in a single interferometer but above pmax = 13 - where the two compared prob- 
abilities are about equal to 50% - it becomes more likely to trigger in two detectors. Finally, the probability of full 
coincidence in the three interferometers of the network is considered. It is quite small for pmax < 10 and reaches 
useful levels only for very high maximal signal to noise ratios: the 50% efficiency is reached only for pmax — 30. 

To conclude this section. Figure nfl presents the detection probability in at least one interferometer of the network 
("OR" strategy) as a function of the source location in the sky. The maximal signal to noise ratio is set to 10 and 
so this graph can be compared with the top map of Figure H. The better efficiency achieved (about double, 67% in 
average instead of 34%) indeed corresponds to a more homogeneous detection probability over the sky. 

In this section and the next more details are given on the coincidence probabilities as simultaneous detection in 
different interferometers will ensure much higher confidence levels for candidate events. 
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C. Twofold coincidences 

Figure 111] shows the detection efficiencies for the three combinations of two detectors: VIRGO-Hanford, VIRGO- 
Livingston and Hanford-Livingston. By construction, the association between the two LIGO antennas is always more 
efficient than VIRGO and one of the LIGO interferometers. Nevertheless, the detection efficiency remains small: it is 
only 20% for the two LIGO antennas with Pmax=10. One can therefore conclude that simultaneous detection in two 
given detectors is unlikely for weak GW signals and becomes likely (« 40%) only for pmax > 15. The last graph of 
this figure presents the twofold coincidence detection probability in the network, when triggering occurs in at least 
two interferometers among the three. For the smallest values of Pmax it corresponds to the two LIGO antennas case 
but then it increases more quickly: already at pmax = 10 this probability is about 50% larger than for the two LIGO 
interferometers. 

As the coincidence detection efficiency remains low for small values of Pmax it is instructive to investigate how it 
is distributed over the sky. Figure ^ presents the detection efficiency of twofold coincidences versus (a, sin J) for 
Pmax = 10 - a likely value for the SNR of the first GW events detected. The three first maps are for the different 
twofold coincidences in the network while the fourth shows the detection probability in at least two detectors. One can 
see that despite the differences in the colour code the areas of high coincident detection probability are of small extent 
for the three configurations. In the last graph these regions are connected and the detection efficiency is higher, but 
nevertheless a major part of the sky remains invisible for the two-detector coincidences for low - but unfortunately 
realistic - pmax values. 

D. Threefold coincidences 

From the conclusions of the previous sections it is clear that a simultaneous detection in the three detectors is not 
likely unless the optimal signal to noise ratio is large. Therefore it appears difficult to reconstruct GW astrophysical 
informations from a source at the expected sensitivity level provided by the first generation of detectors. 

Nevertheless, Figure |l^ shows the detection efficiency sky maps for values of the maximal signal to noise ratio of 
10, 15 and 20 respectively. In the ffist case (pmax — 10), non zero efficiency is concentrated in two small regions 
corresponding to the common visible areas of the LIGO- VIRGO, but even in these parts of the sky the detection 
efficiency is lower than 30%. In the two other cases, the distribution is more uniform with unfortunately some large 
areas which remain invisible - though they decrease as Pmax increases. 
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V. TIMING RECONSTRUCTION 

Determination of the absolute timing of a detected GW burst is an important question. Firstly, coincidences in 
different antennas require time correlations between the respective detected signals. Secondly, the reconstruction of 
the source location is only based on arrival times detected in three detectors located in different places on Earth 
and the better the precision, the smaller the angular error box in the sky. Thirdly, coincidences with other types of 
detectors will also be based on time correlations. An application of the later using neutrinos can be found in [ p5[ . 

Let At be the time difference between the actual arrival time of the GW on the detector trcai and the reconstructed 
time idctoct corresponding to a detection with a given filter in the interferometer's output, when the SNR is maximal. 

At = iioal — ^detect (5-1) 

Generally speaking, there could be an offset between trcai and tdotoct specific to the filter algorithm but here we 
assume it to be zero as would be obtained with matched filtering whose output peaks when the signal and the filter 



overlap |34|. Then the only significant parameter is the RMS of the distribution At^'^^. 



A priori At can be split into two parts: 
• Atsampiing duc to thc discrctc data sampling of the experiments. For a signal of characteristic duration r 3> -J- 



it is completely negligible as it is well approximated by an uniform distribution in thc range 



1 . 1 

2/o ' 2/o 



whose 



standard deviation is J-^ « 1.4 x 10 ^ ms e.g. for VIRGO. But for signals of small duration it can become 
significant by randomly dropping high amplitude parts of the GW whose consequence would be to trigger off 
maximum (or to loose the event). 

• Atnoisc^ if the GW signal is embedded in detector noise, the precise location of the output highest value will 
depend on the actual noise time series; one expects this component to dominate for large r. From dimensional 
analysis it is clear that At must be proportional to r and so it is convenient to use the dimensionless quantity 

Af 

T 

Figure |lj shows the normalised RMS ^^ — as a function of Pmax for different values of the Gaussian width r. 
Apart from the case with t = 0.1 ms which is sensitive to the finite sampling frequency as previously mentioned, all 
the other curves are identical. 
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From the data analysis point of view, it is mandatory to link the timing error with a measured quantity, the filter 
output p. The dependence of At^"^^ on this variable for r values insensitive to the sampling rate can be represented 

by 



At^Ms ^ li^ 



VI ms/ 



P 
The validity range of this equation is r ^ 0.2 ms and the filter output p ^ 5 — 6, as shown on Figure ^ for a particular 

case, T ~ 1 ms and pmax = 20: the fit curve and the real one are in good agreement from p = 6. 

For T = 1 ms and a filter maximal output of about 10, one has At'^^^ k, 0.1 ms which is well below the millisecond, 

the minimum level of precision suitable for coincidences with for instance neutrino detectors |36[| . 

VI. CONCLUSIONS 

This paper deals with the basic aspects of burst coincidence detection in a network of three gravitational wave 
antennas: VIRGO, LIGO Hanford and LIGO Livingston. It docs not focus on the signal parameter reconstruction - 
the inverse problem - but rather on studying the coincidence detection probability to see whether this kind of event 
is likely or not. In this respect, the results obtained are somewhat disappointing: the detection probability in a given 
interferometer is strongly reduced by the beam pattern functions - on the average only about 40%. We also show that 
there is a good complementarity between VIRGO and the two LIGO interferometers for detection ( the 'OR' strategy 
efficiency is higher than 50% for an optimal signal to noise ratio pmax ^ 8), but coincidences are very unlikely for 
weak signals. 

Concerning twofold coincidences, the two LIGO antennas configuration has a much better detection efficiency 
compared to VIRGO with one of the LIGO interferometers (a factor of two larger for pmax = 10). Nevertheless, adding 
VIRGO and looking for any twofold coincidences allows the region of likely detection to be significally extended and 
the efficiency to be correspondingly increased: a 50% enhancement is obtained for pmax = 10, leading to a value 
averaged over the sky slightly higher than 30%). Finally, threefold coincidences between detectors are quite unlikely 
below the value Pmax = 30 where the 50% efficiency is reached. 

As far as timing is concerned, the situation is quite satisfactory: Monte-Carlo simulations show that the RMS 
timing can remain below 1 ms even for low values of Pmax, which is an interesting point for coincidence with other 
kind of detectors. Even if the detection itself remains unlikely, the timing will be accurate enough if the GW is seen 
by the detectors. 
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The basic aspects covered in this paper about beam pattern functions and timing accuracy apply equally well to 
all transient sources whose characteristic times are much shorter than a day; in particular, this is the case for the 
coalescing binary signals which are expected to last from one second to a few minutes in the detector bandwidths. 
So coincidence analysis of binary inspirals will face the same detection efficiency problem as discussed for the bursts. 
Apart from this inescapable fact the situation of the data analysis is different here since the known binary signals 
can be searched for through matched filtering. This opens the possibility of "coherent" rather than "coincidence" 
analyses, as proposed in |l9,M, which are shown pl| to give better results than a simple coincidence analysis; thus, 
this method may compensate for part of the geometric effects presented in this paper by using in an optimal way all the 
information available in the detectors. This improvement is made at the cost of a large increase in computing power 
since the source sky coordinates have to be included in the template definition. Therefore, despite the fact that it is 
sub-optimal, coincidence analysis may be the only available tool in the first years of operation of the interferometers, 
even for binaries and especially for poorly modelled sources such as GW bursts. 

One can a priori imagine two main strategies about their search: 

• if the detector understanding is high enough to allow a proper elimination of non-stationary noise events, a burst 
detection is quite likely in at least one interferometer of the network since they cover the sky in a complementary 
way. Extending the number of sensitive antennas would further improve the detection efficiency. This situation 
is nevertheless not likely to arise in the first periods of data taking; it could even then be hopeless if the 
interferometer noise is not Gaussian. However, this strategy is well-adapted to coincidences with non-GW 
detectors if they have a negligible false alarm rate |Q ; 

• GW detection will in general require coincidences between interferometers. They may not be very frequent 
unless improvements are performed in the detector noise levels. Nevertheless specific strategies must be devised 
in order to maximise sky coverage and detection efficiency p7| . 

The results presented in this article give an overview of the GW burst detectability in a network of antennas 
focusing only on the detection efficiency. To go beyond, it is necessary to take into account the corresponding false 
alarm rates and to compare all the available associations of detectors sensitive to GW or other types of radiation 
such as neutrinos. Forthcoming papers will deal with these questions and also present new ideas about the "inverse 
problem" [p6p|. 
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FIG. 1. Comparison of the averaged beam pattern maps averaged over the polarisation ^ for the three detectors: VIRGO, 
LIGO Hanford and LIGO Livingston. 
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FIG. 11. Twofold coincidence detection probability. Bottom plot: full graph with pmax £ [5; 50]; top plot: zoom on small 
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